Quantum Cloning by Cellular Automata 
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We introduce a quantum cellular automaton that achieves approximate phase- covariant cloning 
of qubits. The automaton is optimized for 1 — »• 2N economical cloning. The use of the automaton 
for cloning allows to exploit different foliations for improving the performance with given resources. 



Quantum cellular automata (QCAs) have attracted 
considerable interest in recent years pQ , due to their ver- 
satility in tackling several problems in quantum physics. 
Quantum automata describing single particles corre- 
spond to the so-called quantum random walks [2] , whose 
probability distributions can be simulated in with an op- 
tical setup H2IU- Solid-state and atom-optics systems, 
such as spin-chains, optical lattices, or ion chains, can be 
viewed as implementations of QCAs, though in a Hamil- 
tonian description. Recently QCAs have also been con- 
sidered as a model of quantum field theory at the Planck 
scale [5j[6]. 

In this scenario, general coordinate transformations 
correspond to foliations, such as those introduced in Ref. 
[7] for operational structures, e.g. the digital equivalent 
of the relativistic boost is given by a uniform foliation 
over the automaton [8] . Besides the link with fundamen- 
tal research, the possibility of foliations makes the QCA 
particularly interesting also for implementing quantum 
information tasks. In this Letter we will explore such a 
potentiality for the case of quantum cloning as a sam- 
ple protocol. We will show how economical cloning can 
be implemented by a quantum automaton, and how the 
foliations can be optimized and exploited for improving 
the efficiency of the protocol. 

It is well known that quantum cloning of nonorthogo- 
nal states violates unitarity [9] or linearity [10 of quan- 
tum theory, and it is equivalent to the impossibility 
of measuring the wave- function of a single system [TT] . 
However, one can achieve quantum cloning approxi- 
mately, for a given prior distribution over input quantum 
states. For uniform Haar distribution of pure states the 
optimal protocol has been derived in Ref. [T2] , whereas 
for equatorial states it has been given in Refs. [T3l fT4] . 
In the present Letter we consider specifically this sec- 
ond protocol, corresponding to clone the two-dimensional 
equatorial states 

|0) = -L(|O)+e^|l)). (1) 

The cloning is phase- covariant in the sense that its per- 
formance is independent of </>, i.e. the fidelity is the same 
for all states \<j>). For certain numbers of input and out- 
put copies it was shown that the optimal fidelity can be 
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Figure 1. Realisation of one-dimensional quantum cellular 
automaton (QCA) with a structure composed of two layers 
of gates A and B. This is the most general one-dimensional 
automaton with next-nearest neighbour interacting minimal 
cells. 

achieved by a transformation acting only on the input 
and blank qubits, without extra ancillas [HI [16]. Since 
these transformations act only on the minimal number of 
qubits they are called "economical" . The unitary opera- 
tion U pcc realising the optimal 1—^2 economical phase- 
covariant cloning is given by [15] 

U pcc \0)\0) = |0)|0), (2) 

u pcc \i)\o) = -L(|o>|i) + |i)|o», 

where the first qubit is the one we want to clone, while 
the second is the blank qubit initialised to input state |0). 
In Ref. [16] the economical map performing the optimal 
N — » M phase-covariant cloning for equatorial states of 
dimension d is explicitly derived for M = kd + N with 
integer k. 

We will now analyse a QCA implementation of the 
economical quantum cloning. Our automaton is one- 
dimensional, and a single time-step corresponds to a uni- 
tary shift-invariant transformation achieved by two ar- 
rays of identical two-qubit gates in the two-layer Margo- 
lus scheme [T reported in Fig. [I] Due to the locality of 
interactions, information about a qubit cannot be trans- 
mitted faster than two-systems per time step, and this 
corresponds to the cell (qubit) "light-cone" made of cells 
that are causally connected to the first. Every event out- 
side the cone has no chance to be influenced by what 
happened in the first cell, thus the quantum computa- 
tion of the evolution of localized qubits is finite for finite 
numbers of time-steps. 

We now remind the concept of foliation on the gate- 
structure of the QCA [8]. Usually in a quantum circuit- 
drawn from the bottom to the top as the direction of 
input output-one considers all gates with the same hor- 
izontal coordinate as simultaneous transformations. A 
foliation on the circuit corresponds to streching the wires 
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Figure 2. (Colors online) Foliations over the automaton. Two 
leafs of two different uniform foliations are depicted with dash- 
lines in different colors (the complete foliation is obtained 
upon repeating vertically the leaf). The systems along each 
leaf are taken as simultaneous. The red "cut" is usually re- 
ferred to as the rest-frame foliation. 



(namely without changing the connections), and consid- 
ering as simultaneous all the gates that lie on the same 
horizontal line after the stretching. Such horizontal line 
can be regarded as a leaf of the foliation on the circuit 
before the stretching transformation. Therefore, a foli- 
ation corresponds to a specific choice of simultaneity of 
transformations (the "events"), namely it represents an 
observer or a reference frame. Examples of different folia- 
tions are given in Fig. [2] Upon considering the quantum 
state at a specific leaf as the state at a given time (at the 
output of simultaneous gates), different foliations corre- 
spond to different state evolutions achieved with the same 
circuit. Therefore, in practice we can achieve a specific 
state belonging to one of the different evolutions, by sim- 
ply cutting the circuit along a leaf, and tapping the quan- 
tum state from the resulting output wires (the operation 
of "stretching" wires should be achieved by remembering 
that by convention the wires represent identical evolu- 
tions, not "free" evolutions). 

We will now show how to perform a 1 4 2N phase- 
covariant cloning of the equatorial states ([I]) with a QCA 
of N layers, with all gates identical, performing the uni- 
tary transformation denoted by A, acting on two qubits. 
Due to causality, we can restrict our treatment to the 
light-cone centered in the state to be cloned and ini- 
tialise all blank qubits to |0), as shown in Fig. [3] By re- 
quiring phase-covariance for the cloning transformation, 
the unitary operator A must commute with every trans- 
formation of the form P x 0P X , where P x is the general 
phase-shift operator P x = exp[|(l — o- z )x] f° r a single 
qubit, with a z the Pauli matrix along z. Therefore, we 
impose the condition 



[A,P x 0P x ] =0, V X 



(3) 



This implies that the matrix A must be of the form A = 
diag(l, V, 1), where V is a 2 x 2 unitary matrix. Notice 
that the transformation A then acts non-trivially only on 
the subspace spanned by the two states {|01), |10)} and 
it is completely specified by fixing V. 

In order to derive the optimal cloning transformation 
based on this kind of QCA we will now maximise the av- 



Figure 3. Cone of gates which contribute to the phase- 
covariant cloning, given the input state \<fi) at site N. 



erage single-site fidelity of the 2Af-qubits output state 
with respect to the unitary operator A. In order to 
achieve this, we write the initial state of 27V qubits in 
the following compact form 



l*0> 
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(4) 



where we define \Q) = |0...0) as the "vacuum state" 
with all qubits in the state, and \k) = |0 ... 01^0 ... 0) 
as the state with the qubit up in the position fc, and all 
other qubits in the state down. Without loss of generality 
in the above notation the qubit to be cloned is supposed 
to be placed at position N and it is initially in the state 
Since the gate A preserves the number of qubits 
up [18], the evolved state through each layer will belong 
to the Hilbert space spanned by the vacuum state and 
the 2N states with one qubit up. The whole dynamics of 
the QCA can then be fully described in a Hilbert space of 
dimension 2N+1. The output state can thus be generally 
written as 



2N 



\* 2N ) = -=(|Q) 



(5) 



k=l 



where the amplitudes of the excited states depend 
only on the explicit form of the gate A. 

The reduced density matrix pk of the qubit at site k 
can then be derived from the output state ([5| as 



Pk 



--Tr- k [\y 2N )(y 2N \ 



1 

2 L 



(6) 
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where Tr^ denotes the trace on all qubits except qubit k. 
The local fidelity of the qubit at site k with respect to 
the input state \<p) then takes the simple form 



Fh = 



\Pk\ 



= -(l + ReK}). (7) 



As we can see, F k depends only on the amplitude of 
the state with a single qubit-up exactly at k. Since the 
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Figure 4. (Colors online). Chromatic map of the local fideli- 
ties in terms of the considered qubit and the layer. The orange 
colour is brighter for increasing local fidelity. The simulation 
involves a number of layers N = 40, while the total number 
of qubits is 2 iV, since it doubles at each layer. 



gate A is generally not invariant under exchange of the 
two qubits, the fidelities at different sites will be in gen- 
eral different. We will then consider the average fidelity 



2N 



J2k=i-Fk as figure of merit to evaluate the per- 
formance of the phase-covariant cloning implemented by 
QCA. Notice that the whole procedure corresponds to a 
unitary transformation on the 27V-qubit system, without 
introducing auxiliary systems, namely it is an economical 
cloning transformation. 

The calculation of the amplitudes was performed 
numerically by updating at each layer the coefficients of 
the state Notice that the amplitude of layer j and 
site k influences only the amplitudes of the subsequent 
layer j + 1 and sites either k — 1, fc or fc, A; + 1, depending 
on whether the state |1) enters in the right or left wire of 
A. The action of A on the qubits j and j + 1 is given by 



A(jJ + l)\k)={ 



V22 1 j) +v 12 \j + 1) if k = j 
V2i\j)+vn\j + l) iffc = j + l (8) 
\k) otherwise , 



where Vij are the entries of the operator V in the basis 
{|01), 1 10)}. Notice that the vacuum state \Q) is invariant 
under the action of A. The iteration of Eq. for each 
layer then leads to the amplitudes of the output state ([5| . 

As a first explicit example we will consider a QCA 
employing the optimal 1 — >• 2 phase-covariant cloning 
([2]). In this case the gate A must implement the unitary 
transformation ([2| . The non-trivial part V of gate A can 
then be chosen to be 
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(9) 



where all coefficients are real. The corresponding lo- 
cal fidelities at every layer are reported in Fig. [4] As 
we can see, the figure exhibit fringes of light and dark 
colour. Moreover, the light-cone defined by causality can 
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Figure 5. (Colors online) The average 1 — >> 2N phase- 
covariant cloning fidelity achieved by the QCA in the rest 
frame (see Fig. [2|. Comparison with the optimal econom- 
ical phase-covariant cloning in Ref. [14 . The purple dots 
represents the QCA cloning optimized over the gate unitary 
A. The blue dots correspond to the use of gate A achieving 
the optimal 1—^2 cloning U pcc in Eq. |2]). The yellow dots 
represent the unrestricted optimal economical cloning. 



be clearly seen: outside this cone no information about 
the initial state can arrive, thus every system has the 
same fidelity of 1/2. Notice that there is a sort of line, 
approaching the right top corner, along which the fidelity 
is quite high. Regarding the local fidelities of the final 
states, they are in general quite different from each other, 
and can vary very quickly even between two neighbour- 
ing qubits. The average fidelity is reported in Fig. [5] as 
a function of the number of layers. Notice that the av- 
erage fidelity of the optimal economical phase-covariant 
quantum cloning (without the constraint of automaton 
structure) approaches the value 3/4 for a large number 
of output copies [T4] . 

In order to improve the average fidelity we have then 
maximised it with respect to the four parameters defin- 
ing the unitary operator V. Numerical results achieved 
up to N = 20 show that the optimal cloning performed 
in this case is not much better than the one given by 
the iteration of The latter turns out to be outper- 
formed only when the number of layers composing the 
automaton is odd, as shown in Fig. [5] Further numerical 
results show that no gain can be achieved if the automa- 
ton is composed of layers of two different gates A and B. 
Actually, in this case it surprisingly turns out that the 
optimal choice corresponds to B = A, namely we do not 
exceed the average fidelity obtained by employing a single 
type of gate. Since all one-dimensional QCA with next- 
nearest neighbour interacting cells with two qubits can 
be implemented by a two-layered structure, we have then 
derived the optimal phase-covariant cloning transforma- 
tion achievable by the minimal one-dimensional QCA. 

We will now show that the average fidelity can be im- 
proved by considering different foliations. Suppose that 
we are given a fixed number M of identical gates A to im- 



4 




Figure 6. Illustration of the classification of foliations. A 
possible foliation with M — 6 gates is given. From the corre- 
spondence between the gates lying under the "cut" and the 
rotated dots on the right, we identify this foliation with the 
partition {4, 1, 1}. 



M(Layers) 


Frest F 


Optimal Foliation 


1(1) 


0.853 0.853 


{1} 


3(2) 


0.676 0.693 


{3} 


6(3) 


0.617 0.679 


{2,2,2} 


10(4) 


0.588 0.670 


{4,3,3} 


15(5) 


0.570 0.653 


{4,4,4,3} 


21(6) 


0.558 0.614 


{4,3,2,2,2,2,2,2,2} 


28(7) 


0.550 0.603 


{6,6,6,5,5} 



Table I. Results of the maximisation over foliations up to 
M — 28, i.e. QCA composed of up to 7 layers. 

plement a QCA. We are then allowed to place the gates 
in any way such that the causal structure of the consid- 
ered automaton is not violated. Which is the configura- 
tion, i.e. the foliation, that performs the optimal phase- 
covariant cloning for fixed Ml In this framework we have 
to maximise not only over the parameters that define V 
but also over all possible foliations. Thus, the M fixed 
gates play the role of computational resources, and the 
optimality is then defined in terms of both the parame- 
ters characterising the single gate A and the disposition 
of the gates in the network. As a first example, suppose 
that we are given M = 3 gates. In this case there are 3 in- 
equivalent foliations: one for the rest frame (see Fig. [2]), 
and two along the straight lines defining the light-cone. 
As expected, for increasing M the counting of foliations 
becomes more complicated and the problem is how to 
choose and efficiently investigate each possible foliation. 
It turns out that the problem of identifying all possible 



foliations of a QCA of the form illustrated in Fig. [3] for a 
fixed number of gates M is related to the partitions of the 
integer number M itself (by partition we mean a way of 
writing M as a sum of positive integers, a well known con- 
cept in number theory [19]). Two sums that differ only 
in the order of their addends are considered to be the 
same partition. For instance the partitions of M = 3 are 
exactly 3 and given by {3}, {2,1}, and {1,1,1}, while the 
partitions of M = 6, corresponding to a 3-layer setting 
in the rest frame, are 11 and given by {6}, {5, 1}, {4, 2}, 
{4,1,1}, {3,3}, {3,2,1}, {3,1,1,1}, {2,2,2},{2,2,1,1}, 
{2, 1, 1, 1, 1}, and {1, 1, 1, 1, 1, 1}. The link between folia- 
tions and partitions is illustrated in Fig. [6j which shows 
how partitions can be exploited to identify foliations. For 
a fixed number of gates M the number of foliations is then 
automatically fixed and each foliation corresponds to a 
single partition. The correspondence is obtained as fol- 
lows. Each addend represents the number of gates along 
parallel diagonal lines, starting from the vertex of the 
light-cone, as shown in Fig. [6] for the particular case of 
M = 6. 

Based on this correspondence, we can investigate the 
performance of the phase- covariant cloning as follows. 
For a fixed foliation, we first maximise the average fidelity 
with respect to the gate A. Then we choose the highest 
average fidelity by varying the foliation. We worked out 
this procedure numerically for M = 1, 3, 6, 10, 15, 21, 28, 
i.e. the number of gates composing the QCA with 
N = 1,2,3,4,5,6,7 layers, respectively. Our results are 
shown in Table [T| where the maximisation in the rest 
frame is also reported for comparison. As we can see, 
exploiting different foliations leads to a substantial im- 
provement of the average fidelity. 

In summary, we have introduced a new way of achiev- 
ing quantum cloning through QCAs. We have derived 
the optimal automaton achieving economical phase co- 
variant cloning for qubits. We have shown how the fi- 
delity of cloning can be improved by varying the foliation 
over the QCA, with fixed total number of gates used. By 
developing an efficient method to identify and classify fo- 
liations by means of number theory, we have optimised 
the performance of the QCA phase covariant cloning for 
a given fixed number of gates, and obtained in this way 
the most efficient foliation. 
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